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ABSTRACT 

Pulsars and central engines of long gamma ray burst - collapsars - may 
produce highly magnetized (Poynting flux dominated) outflows expanding in a 
dense surrounding (interstellar medium or stellar material). For certain injection 
conditions, the magnetic flux of the wind cannot be accommodated within the 
cavity. In this case, ideal (non-dissipative) MHD models, similar to the Kennel 
and Coroniti (1984) model of the Crab nebular, break down (the so-called a 
problem). This is typically taken to imply that the wind should become particle- 
dominated on scales much smaller than the size of the cavity. The wind is then 
slowed down by a fluid-type (low magnetization) reverse shock. Recent Fermi 
results, indicating that synchrotron spectrum of the Crab nebula extends well be- 
yond the upper limit of the most efficient radiation reaction-limited acceleration, 
contradict the presence of a low sigma reverse shock. 

We propose an alternative possibility, that the excessive magnetic flux is de- 
stroyed in a reconnection-like process in two regions: near the rotational axis and 
near the equator. We construct an example of such highly magnetized wind hav- 
ing two distinct reconnection regions and suggest that these reconnection cites are 
observed as tori and jets in pulsar wind nebulae. The model reproduces, qualita- 
tively, the observed morphology of the Crab nebula. In parts of the nebular the 
dissipation occurs in a relativistically moving wind, alleviating the requirements 
on the acceleration rate. 



1. Introduction 

The interaction of strongly magnetized, relativistic outflows with a dense surrounding 
is a basic problem in pulsar physics, which recently became important for gamma ray bursts 
(GRBs) central engines. This seemingly straightforward problem turns out to be a compli- 
cated one, as is illustrated by what became known as the "a problem" (Kennel & Coroniti 
1984). Following the work of Rees & Gunn (1974), Kennel & Coroniti (1984) assumed that 
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the central source produces a wind which blows a cavity in the surrounding medium. The 
key assumption of the model is that in the bulk of the cavity the wind is supersonic, causally 
disconnected from the source. Kennel & Coroniti (1984) found that for a non-relativistic ex- 
pansion of the bubble, i> ou t *C c, if the Poynting flux in the wind is comparable to or exceeds 
the particle flux, there is no acceptable solution for the structure of a resulting cavity. Intro- 
ducing a parameter a as the ratio of Poynting -F Poynting to particle F p fluxes, a = F Poynting / F p , 
they found that solutions where wind is injected into cavity with relativistic velocities, while 
the cavity expands non-relativistically, exist only for a ~ v out /c <C 1. Formally, for a > 1 
the reverse shock in the wind is relativistic with respect to the downstream flow, and since 
the flow matches to the slowly expanding cavity, the shock is driven back into the pulsar on 
a light crossing time. As a result, the initial set up of the problem, that of a central engine 
producing relativistic wind, breaks down. The requirement a <C 1 in the bulk of the cavity 
runs contrary to models of pulsar magnetospheres, which all predict generation of highly 
magnetized, a ^> 1 wind (Goldreich & Julian 1970). These expectations are confirmed by 
recent numerical simulations (Spitkovsky 2006). Kennel & Coroniti (1984) concluded that a 
must change from being much larger than unity at the light cylinder to much smaller than 
unity in the bulk of a PWN. 

The conversion of high a to low is very hard to achieve in ideal MHD (Heyvaerts & 
Norman 2003), though a number of recipes have been proposed (e.g. Vlahakis 2004). It 
is beyond the scope of this paper to discuss this mathematically very difficult problem, but 
simple, non-dissipative models do not achieve high to low a conversion (the applicability of 
fluid approach has also been questioned Lyubarsky 2003). 

Note, that the requirement of evolution of a from high to low follows from a self- 
consistency of the model, in a sense that for a > 1 the whole model, that of a central pulsar 
producing a supersonic wind which is shocked within a PWN, becomes inconsistent. Thus, 
a question still remains, if a remains high, how would a flow behave?. This is the main 
question addressed in this paper. 

Recently this problem became important for long GRBs. Detection of Type Ic super- 
novae nearly coincidence with long GRBs unambiguously linked them with deaths of massive 
stars (Stanek et al. 2003; Hjorth et al. 2003). Studies of the host galaxies of long GRBs, 
which turned out to be actively star-forming, further strengthens this association (Djorgov- 
ski 2001). The outflow resulting in a GRB can be driven by a "millisecond magnetar" (a 
fast rotating strongly magnetized protoneutron star, e.g. Usov 1992) or a black holes with 
an accretion disc (e.g. MacFadyen & Woosley 1999). In both cases the wind produced by 
the central engine may be Poynting flux dominated, u»L The evolution of the resulting 
bubble inside the star then resembles the evolution of a PWN. Again, in order to achieve 
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a self-consistent model, it is required that a changes from high to low values (Barkov & 
Komissarov 2008; Bucciantini et al. 2009). (In passing we note that the electro-magnetic 
model of GRBs (Lyutikov & Blandford 2003; Lyutikov 2006) explores the possibility that a 
remains high after the flow breaks out from the star.) 



2. Radiation reaction-limited acceleration of leptons and Fermi spectrum of 

Crab 

There is an upper limit on the frequency of synchrotron emission by radiation reaction- 
limited acceleration of electrons. If the accelerating electric field is a fraction rj < 1 of the 
magnetic field (this is equivalent to acceleration on time scale of inverse cyclotron frequency 
l/(r]UB,rei), where u>B,rei = 7/Vb is relativistic cyclotron frequency of a particle), equating 
the acceleration rate and synchrotron energy losses, 

rjeBc^—^Ws (1) 
the peak frequency of emission becomes 

7TI C 

lo = r)u max , u max = 0.44^ r (2) 

This corresponds to energies (see also de Jager et al. 1996). 

mc 

e max = Huj max = 0.44ft— « 30 MeV. (3) 

e l 

Note, the upper limit (3) assumes non-stochastic, DC-type acceleration. The break energy 
observed by Fermi satellite in Crab is ~ 100 MeV (Fermi Collaboration & Pulsar Timing 
Consortium 2009), already three times higher than the limit (3). If the spectrum above 
the break requires presence of higher energy particles, those particles should be accelerated 
by an electric field much larger than the magnetic field. This is highly problematic under 
astrophysical conditions, where abundant supply of particles ensures that E < B. 

There is a number of possible resolutions of this contradiction. First, magnetic field can 
be strongly inhomogeneous, spatially or temporarily, so that particles are accelerated in a 
low magnetic field regions, e.g. very close to the equatorial current sheet, and then emit 
when they enter a high field region. The possibility of spatially inhomogeneous magnetic 
field is also consistent with the results of recent numerical simulations of relativistic shocks, 
indicating that only unmagnetized shocks are efficient accelerators (Sironi & Spitkovsky 
2009, another case of highly efficient acceleration, quasi-parallel shocks, is not applicable to 
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PWNe). The problem with this scenario is that the Crab nebular is a very efficient emitter: 
about 20% of the spin-down luminosity is converted into radiation. If acceleration is limited 
to a narrow current layer containing only a small part of the spin-down luminosity, this put 
an unreasonable demands on the efficiency of acceleration. 

Temporal variations in the strength of magnetic field near the termination shock, seen 
in numerical simulations by Camus et al. (2009), remain is distinct possibility, with a clear 
prediction: the Crab synchrotron spectrum should be strongly variable on time scale of 
month, of the order of wisp variability, with the cut-off energy dropping well below 100 MeV. 
This can be tested with the ongoing observations by Fermi satellite. Previously, EGRET 
data did indicate a moderate level of the cut-off energy variability (de Jager et al. 1996). 

Another possibility to explain a high acceleration rate, is that the flow is relativistic in 
the post-shock region. In this case the observed emission is beamed by the bulk Lorentz 
factor r, oj ~ r]Tuj majX . Thus, a high post-shock Lorentz factors T may alleviate the problem 
of super-efficient acceleration. Relativistic post-shock flows are indeed expected within the 
Komissarov & Lyubarsky (2004) model of PWNs (see also Del Zanna et al. 2004). In 
this model, the reverse shock is generally oblique, so that the post-shock flow is moving 
relativistically. But the post-shock flow is only mildly relativistic, independently of the pre- 
shock Lorentz factor, r ~ 1/6 ~ few, where 9 ~ 1, is the angle that the radially moving 
unshocked wind makes with the forward shock (this angle can be small near the axis, but 
very little energy is coming along those directions). 

As we show below, the high sigma flows allow dissipation in a relativistically moving 
plasma, boosting the maximum observed synchrotron frequency to w ~ r]Tu max and allevi- 
ating the requirements on r], which can be of the order r] ~ 1/r <C 1. 

3. The a problem 

The a problem of Kennel & Coroniti (1984) implies that ideal, non-relativistic, homol- 
ogous expansion of a bubble of strongly magnetized plasma injected with nearly a speed of 
light, cannot occur. The reason is that magnetic flux and energy are supplied to the inflat- 
ing bubble by the rates that cannot be accommodated in the bubble. The rate of supply is 
determined by the processes inside the light cylinder of a compact object, while inflation of 
the bubble is controlled by the external gas density. Formally, the surface of the bubble will 
influence back on the source and would start to control the injection rate, consistent with 
the allowed rate of expansion. But even if the wind remains subsonic, it is unlikely that 
processes at the edge of the inflating bubble would influence the wind generation region near 
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the light cylinder due to large separation in scales, thousands in case of GRBs and billions 
in case of PWNs. 

One way to approach the a problem, is to try to solve magnetohydrodynamical reverse 
shock conditions in the strongly magnetized wind, that match the slowly expanding PWN 
cavity downstream. This can only be done for low er: for a > 1 the reverse shock is driven 
back into the pulsar on light crossing time at which point the model becomes meaningless 
(Kennel & Coroniti 1984). 

To illustrate the a problem from another viewpoint, assume that the central source 
injects a relativistic wind that carries electromagnetic power and magnetic flux in a form of 
toroidal magnetic field. The central source injects energy and magnetic flux into a spherical 
cavity, PWN, with the injection velocity Vi n = f3i n c and the outer boundary of the cavity 
expands with velocity v ou t = cf3 ou t <C v « n . The core assumption here is that the injection ve- 
locity is independent of the conditions in the nebular. We assume that the expansion velocity 
of the nebular is much smaller than both the injection velocity and the fast magnetosonic 
speed inside a bubble, so that the interior is in a state of a quasistatic equilibrium. 

Let us assume that magnetic flux is injected into the Crab nebular at some rate $. 
This rate depends on the integrated polar angle-dependent luminosity L{9) produced by 
the pulsar. Inside the nebular, where the velocity of expansion becomes strongly sub- 
fastmagnetosonic, the magnetic field will relax to a minimum energy state consistent with 
a given magnetic flux. It may be shown that such a state corresponds to the current con- 
centrated on the symmetry axis. Since such a current distribution would correspond to the 
infinite current density, we will assume that there is a current-carrying core of conical shape 
with opening angle 9 C <^ 1. Equivalently, the cavity is in a state of dynamical balance and 
since the pressure forces are negligible, the magnetic field relaxes to a state of (quasi)static 
force-free equilibrium with no velocity shear and associated electric fields. The only possible 
equilibrium of a purely toroidal magnetic field corresponds to a line current (see, e.g. , Eq. 
(13)), with infinite current density on the axis. This state should have the magnetic flux 
equal the total injected flux. 

The central source may have anisotropic energy flux. As an exemplary case, consider 
power distribution produced by an aligned pulsar Michel (1973). In this case L oc sin 2 9, and 
normalizing magnetic field in the wind to the total luminosity L = 2nrf n (3 in c J j^sin9d9, 
the magnetic field at the inner boundary of the nebular is 




(4) 



The rate of the injection of the magnetic flux is $ = (3 in c J Brdrd9 = 2y/37rP in Lc, so that 
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total flux stored in a cavity is 

$ tot = y/6/3 in Lct (5) 
(and total energy E = Lt). We assume for simplicity that the luminosity is constant in time. 

Inside the nebular, magnetic field corresponds to a current-carrying core of conical shape 
with opening angle 9 C . Equating the total stored flux in such a configuration to the injected 
flux, the magnetic field inside a PWN is 

d y/6cf3 in Lt ( 

in r sin 0i2(t) In tan(0 c /2) w 

where R(t) is the outer radius of a PWN. The ratio of the stored to the injected energy is 

fBl/(8n)dV _3P m 1 

Lt 2/3 ottt lntan# c /2' 1 ' 

where we assumed R(t) = f3 ou tct. For injection velocity of the order of the speed of light, 
Pin ~ 1 and expansion velocity several thousand km/sec, (3 out ~ 1CT 2 , the angle 9 C required 
to accommodate all the injected flux and energy is unphysical small. Thus, the assumption 
of the model, of a source injecting at a given rate magnetically dominated, ideal flow into 
slowly expanding cavity, lead to unphysical consequence. 

Finally, in the simplest example of the paradox, consider a cavity of fixed volume, thus 
neglecting the expansion completely, The energy and the magnetic flux are injected into the 
cavity by the wind linearly in time. But since toroidal flux grows linearly in time, the stored 
magnetic energy grows quadratically, in clear contradiction to linear injection of energy (Rees 
& Gunn 1974). 



3.1. Possible resolutions of cr-paradox 

There are several possible resolutions of this "paradox". First, our assumption of high 
a in the bulk of PWN can be wrong. This is the resolution chosen by Kennel & Coroniti 
(1984), who postulated that that a in the wind decreases from a 3> 1 to a <C 1. 

What if a remains high in the bulk of the flow? One possibility is that the flow would 
break out of the cavity: in laboratory devices, where the sizes of the magnetic field generator 
and the enclosing cavity are similar, and where magnetic field is limited by strength of cavity's 
material, this sometimes leads to spectacular accidents. In case of GRBs this will occur at 
the time of jet break-out, but cannot occur in PWNe. 

Another possibility is that the cavity affects the central engine, so that the rate of 
injection of the magnetic flux adjusts to the rate with which it can be accommodated. We 
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consider the possibility unlikely due to large separation of scales. It is unimaginable that 
a PWN boundary, located at distances of parsecs away from the central source, affects 
generation of pulsar winds at scales of hundreds of kilometers. In the case of GRB central 
engines, the separation of scale is somewhat smaller, only three-four orders of magnitude, 
from 10 7 cm (typical size of light cylinder radius of millisecond neutron star or Schwarzschild 
radius of a black hole) to 10 11 — 10 12 sm (size of a precollapse type Ib/c progenitor), but still 
there are some three-four orders of magnitude. We assume that injection is independent of 
the expansion and investigate consequences of this assumption. 

Another possible resolution of the u-paradox requires destruction of magnetic flux, which 
in turn requires resistivity and reconnection (reconnection in its original sense, as a process 
destroying magnetic flux, not necessarily the magnetic energy). Thus, dissipation must 
become important in the flow. Typically, reconnection is most efficient in a limited regions 
of space associated with high electric current concentration, current sheets. Current sheets 
produced locally, through turbulent motion of plasma, won't destroy the large scale magnetic 
flux. The destruction of the the large scale toroidal magnetic flux can be achieved in two 
generic locations: near the axis in an O-type reconnection, or near the equatorial plane, 
where magnetic field lines of opposite polarity approach each other (assuming, for simplicity, 
an aligned rotator). 

Thus, a flow may be ideal in the bulk and strongly resistive in regions of small measure. 
What would a structure of the flow be in this case? As the magnetic flux and necessarily 
the magnetic energy (and, thus, magnetic pressure) are dissipated, this would create a force 
misbalance that would push the flow towards the dissipation regions. Thus, a flow structure 
brining magnetic field toward dissipation regions - rotational axis and equator - is set up. 
This, in turn, leads to the pile-up of magnetic field near the dissipation regions and to faster 
radial expansion in the direction of those regions (akin to toothpaste tube effect Lyutikov 
& Blandford 2003). This effect of magnetic pile-up would be most pronounced near the axis 
and may be responsible, as we propose, for generation of GRB and pulsar jets. 

Previously, Kirk & Skjaeraasen (2003) discussed magnetic field dissipation in the pulsar 
wind. In their model, the reconnection occurs locally, between different magnetic polarities 
in a striped wind. In this case, the total injected magnetic flux is zero if averaged over the 
pulsar's period. In contrast, we assume that the pulsar rotation axis is nearly aligned with 
the magnetic axis, so that the magnetic field has a given polarity in each hemispheres. 



-8 - 



4. Magnetized flow within a spherical cavity 

Next we consider an idealized problem of a high-a dissipative flow bounded by a cavity. 
We expect that a flow pattern is set up in which a magnetized wind is injected at small 
radii and is deflected at intermediate regions towards the dissipation regions. First, instead 
of calculating the form of the cavity self-consistently, we assume it to be fixed and then 
calculated a flow within it. As the simplest example, we consider a spherical cavity, but the 
following derivation may be repeated for any of the separable coordinates. Second, we con- 
sider a case of highly magnetized plasma a ^> 1 and completely neglect inertial contribution, 
the so-called force-free approximation (Gruzinov 1999, note, that the relativistic force-free 
approximation is different from vacuum, since massless charges create electric currents and 
charge density that affect the flow dynamics). Third, we assume that dissipation is limited 
to an infinitely small volume (formally of measure zero) concentrated near the axis and the 
magnetic equator, while in the bulk the flow is ideal. Fourth, we assume that the flow is 
generated within a volume small compared to the size of the cavity, so that the magnetic 
field is dominated by axially symmetric toroidal field. Neglect of the poloidal magnetic field 
is, perhaps, mathematically the most important simplifying assumption, since it untangles 
the structure of the poloidal current from the structure of magnetic flux surfaces, the main 
complication in solving the Grad-Shafranov equation. Finally, we assume that the flow is in 
a steady state condition, which implies that the velocity of the flow are much larger than 
the velocity of the cavity expansion and that the pulsar is an aligned rotator. 

Equations of ideal force-free electrodynamics (Gruzinov 1999) include Maxwell's equa- 
tions, 

<9E 

— = VxB-47rj (8) 

7* = ~ VxE (9) 

and force-balance 

pE+jxB = 0, (10) 
which allows one to relate the current to the electro-magnetic fields 

. (ExB)V-E+(B-VxB-E-VxE)B , N 

J = 4^ (U) 

In the asymptotic domain, the only nonzero components of the fields are B = B^, E r , E e 
and equations (9) and (11) give 

d t B = --d r (rE e ) + -d e E r 
r r 
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d t E e = --d r (rB) + § (-^deismOEe) + \d r (r 2 E r ] 
r B \ r sin 9 r z 



d t E r = -^—d e {sin6B) - § ( — L-^(sin^) + -^(r 2 ^) ) (12) 
rsmfl B \rsmV r z 

In this formulation, the electrical currents arise exclusively due to charge transport across 
magnetic field by electric drift, j r = p'ff , je = ~P^i where p = V • E/(47r) is charge density. 
Thus, they are not subject to resistive decay in the bulk (e.g. Lyutikov 2003). Toroidal 
axially symmetric magnetic field is = 21 j{cr sin 9), where / is a current flowing through 
a magnetic loop at radius r and polar angle 9. In ideal, time-stationary case the electric field 
is poloidal and irrotational, E = — V$. The force-free condition then gives 

2V/ 2 = r 2 sin 2 #A$V$ (13) 

The neglect of poloidal magnetic flux makes the problem under-determined, in a sense 
that for a given toroidal magnetic field we may chose various distributions of electric field, 
the only constraint being that the flow lines do not intersect the given form of the outer 
boundary. In case of axial symmetry, the force balance gives one scalar equation (13) for two 
functions, toroidal magnetic field and poloidal electric field, or, equivalently, electric current 
and charge density. (Note, that in absence of a shear, $ = 0, Eq. (13) implies / =constant, 
and magnetic field corresponding to a line current.) 

Equation (13) should be solved for a given form of the cavity. If the boundary is 
impenetrable, the velocity of plasma, 

ExB = rsin^xV$) 
B 2 21 1 ' 

should vanish on the boundary. Thus, for a fixed spherical cavity of size R, it is required 
that 

<%$(r = R) = (15) 

Equation (13) is underdetermined since it involves two unknown functions. A partic- 
ularly simple choice of electric potential $ and current I corresponds to vanishing charge 
and current densities in the bulk, so that A$ = and / =constant. (Another analytically 
tractable case is discussed in Appendix A.) Taking the first terms in a series of harmonic 
functions (Legendre polynomials with m — 0, 1), the potential that satisfies boundary con- 
dition (15) is 

/ R\ tan | 
* = *•(!- 7j ' n ta4 (16) 
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where constant 9* is the polar angle that separates flow lines that end on the axis, for 9 < 9* 
and at equator, for 9 > 9*, where 9 is the polar angle at the injection radius r = r Q . The 
flow lines are then given by (see Fig. 1) dgr = —de§/d r &, which integrates to 



In ^Rr 



#ln^# +r ln-^ 



(17) 



The flow lines converge either on the axis or at the equator. It is implicitly assumed 
that there the dissipation sets in and destroys the magnetic field at the end points of flow 
lines. On the axis all the flow lines asymptote to the point 9 = 0, r = R, they never cross 
the line 9 = at r < R. Situations at the equator is very different, a flow line that starts at 
r , 9o reaches equator at 

r R In cot (0*/2) 



tan ^0 

r lncot(fl /2) + Rln-^ 



(18) 



Thus, magnetic energy must be dissipated continuously and effectively at the equator, while 
near the axis a toothpaste-like effect pushes the flow in radial direction. Thus, on the axis 
all the dissipation (formally) happens at a point 9 = 0, r = R. 

The flow velocity corresponding to Eq. (16) is 

v r = A) ( — - 1 
\ r 

n Rsin9 , tan § 

ve = Po In- |r 19 

r tan \ 

where we introduced (3q = $ /(2J). The flow settles down on the axis, vg — > as 9 — > 0, and 
has a finite ^-velocity at equator, v e — > — /5 lntan ff(R/r) at 9 — > ir/2. 

The corresponding Poynting flux toward the equator is 

^ 2(i?-r) 

r r> ■ 2 n 

r d sin y 

„ 2i? , tanf 

In- |rJ$o (20) 



sin^r 3 tan 2 



The energy flux towards the reconnection regions is oc rsm9F g . Near the axis, the energy 
flux diverges as ln#; assuming that resistivity becomes important within a core defined 
by the polar angle 9 C , the ratio of energy fluxes going to the axis and to the equator is 
£ = ln6 l c /lncot(^ /2). We expect that in case of PWNe, the observed relative brightness of 
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Fig. 1. — Flow lines given by eq. (17). The cavity is assumed to be spherical. Injection of 
the wind occurs at some small radius r , labeled NS. Flow lines that star at 9 < 9* (45° in 
this example) terminate near the axis, while those that start at 9 > 9* terminate on equator. 

jet and tori depends on the relative energy fluxes to the corresponding reconnection regions, 
so that the parameter £ is a measure of the relative brightness within this model. At a 
present stage, it is expressed in terms of two formal parameters, the core angle 9 C and the 
separatrix angle 9*. 
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To summarize, we found a flow of strongly magnetized plasma within a given spherical 
cavity. The particular model, chosen mostly for simplicity, has a separatrix angle 9* that 
delineates the flow lines ending on the equator and on the axis. The ratio of energies 
dissipated depends also on the core angle 9 C , where dissipation sets in near the axis. The 
overall value of the velocity scales with the parameter /3 Q . 



4.1. Applicability 

So far, the parameters f3 Q and 9* were arbitrary. We should now check the conditions 
that the resulting solutions correspond to a physically meaningful velocity 

v = Po^(l- R/r) 2 + (R/ry In 2 sin 2 9 < 1 (21) 

at least somewhere inside the sphere r = R. Demanding that v < 1 on the equator at r = it! 
requires that 

n* 

A)lncoty<l (22) 

For parameters satisfying Eq. (22), the condition v < 1 is violated for sufficiently small radii, 

r/R < (3q\J\ + In 2 cot 9*/2. Thus, the model is not applicable at very small radii, where 
mass loading should be important. This is an artifact of our complete neglect of inertia. 



4.2. Surface brightness of a model PWN 

To produce an example of the observed PWN in our model, we assume that (i) the 
emission is generate in the equatorial plane and along the rotation axis; (2) the observed 
intensity is proportional to 5 2 ' 5 , where 5 = 1/ (7(1 — v • n)) is a Doppler factor, 7 is a Lorentz 
factor corresponding to velocity (19), n is a unit vector along the line of sight towards the 
observer; (3) chose inclination angle inferred for Crab nebular, 9 b = 71/ 3 (the same value 
as used by Komissarov & Lyubarsky (2004)); (iv) chose 9* = n/3 (smaller values of 9* 
give large Doppler factors and, correspondingly, more pronounced front-back asymmetry; 
(v) parameter (3 = 0.25 (the values of (3 mostly controls the size of the inner region, where 
the force-free approximation is not applicable; smaller f3 correspond to smaller size of that 
region); (vi) axis inclination of 135° with respect to north-south direction. Results of the 
calculations are presented in Fig. 2. 

Chandra observation of Crab nebular show that the pulsar is located inside a low bright- 
ness bubble, outlined by bright wisps. In our model, the "sudden" onset of dissipation is an 



Fig. 2. — A negative image of the Crab nebular in the high-a model. Point in the center 
denotes the central neutron star. Opening angle of the jets is chosen arbitrarily, 0.1 radian, 
for illustration purposes. Equatorial and jet emission is proportional to local Doppler factor. 
Within the central cavity the model is not applicable; see text for more details. 

effects of Doppler boosting. Dissipation occurs in a limited region near equator and along 
the axis. At small radii, when radial velocity is large, the emission is boosted away from 
the observer. Further out, the flow slows down, as the radial velocity decreases with radius 
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oc R/r — 1, (19) emission from corresponding parts comes into view. The resulting structure 
shows qualitative agreement with observations. 

5. Discussion 

In this paper we outlined a model of a high-a (Poynting flux dominated) flow confined 
by a medium so dense, that the rate of the expansion of the blown-up cavity is not sufficient 
to accommodate the injected toroidal magnetic flux. The model postulates that excessive 
magnetic flux is destroyed within the nebular at a rate dictated by the expansion of the 
nebular. We do not address the related dissipation of energy. Formally, destruction of 
the flux can proceed at an infinitely small resistivity (a fast dynamo model, Vainshtein & 
Zel'dovich 1972, is an example of model with non-conserved magnetic flux, yet negligible 
magnetic energy dissipation). 

The nebular flow in this case is drastically different from the more conventional MHD 
model of PWNe by Rees & Gunn (1974) and Kennel & Coroniti (1984). The assumption of 
the model is that if the injection rate of toroidal magnetic field exceeds the rate at which 
magnetic flux can be accommodated, a flow will find ways to destroy the excessive magnetic 
flux. We present an example of how a nebular can achieve this: a flow pattern is set up, where 
the magnetic energy is injected nearly radially at small distances and is deflected toward 
reconnection cites in the bulk. We give an idealized example of such flows, in which the 
electric current and electric charge density are confined to two regions: near the symmetry 
(rotational) axis and equator. In this particular solution the bulk of the flow is current 
and charge density-free. This solution is very simple and is meant to illustrate a possible 
character of the flow with the nebular. In spite of it simplicity it reproduces qualitatively 
two dissipation regions observed in PWNs, axes and tori, as a well as collimated jet-like flows 
observed in some PWNe. 

The model has implications for dynamics of collapsar winds at early stages of GRB 
explosions. Even in the current simple version of the model, the flow shows a "toothpaste 
tube" effect: the flow is collimated towards and is pushed along the symmetry axis. Our 
assumption of a rigid spherical boundary precluded formation of a jet. Thus, dissipation on 
the axis leads to formation of the flow towards the axis and increased magnetic field. This 
provides an additional (to the ideal case) collimation mechanism. 

In addition, we expect that in the equatorial plane the flow also expands faster: conver- 
sion of the magnetic field energy density into kinetic energy of particles removes the confining 
hoop stresses. Thus we expect that both on the axis and at the equator a PWNs expands 
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faster than at midlatitudes. 

This simple model misses a number of important effects: the properties of the solution at 
small distances do not quite correspond to what is expected from the pulsar injection. In the 
analytical example we discussed the angular dependence of luminosity does not correspond 
to what is expected from pulsar models (e.g. Spitkovsky 2006) (the model has more energy 
near the axis than expected. Also, at small distances the flow is formally superluminal. 
Dissipation of magnetic field inside the reconnection regions should lead to high plasma 
pressure, invalidating the assumption of magnetically dominated plasma. In the current 
formulation, the energy just disappears inside the reconnection regions: this is partly can 
be attributed to radiative losses, but physically most of the energy just change form, from 
magnetic to internal and bulk kinetic, and this will affect dynamics close to reconnection 
regions; the shape of the cavity was assumed to be fixed and given: it should be found 
self-consistently. Finally, we assumed an aligned rotator, while Crab is, in fact, nearly 
orthogonal (Moffett & Hankins 1999). For orthogonal rotators the injected magnetic flux is 
zero if averaged over the period of rotation, so that the reconnection can proceed locally, not 
involving large scales flow discussed in the current model. 

A possible observational distinction between shock acceleration an reconnection is a 
spectral index of accelerating particles p. Shock acceleration typically produces p ~ 2 
(Blandford & Eichler 1987), where dn/d'j oc 7~ p (see, though Lyubarsky & Liverts 2008). 
Non-linear shock acceleration model can, in principle, produce harder spectra, upto p — 1.5 
for adiabatic index of 4/3. Observations of Vela pulsar (Kargaltsev & Pavlov 2004) indicate 
p = 1.1, much flatter than expected. In contrast, acceleration in relativistic reconnection is 
expected to generically produce p — 1 (Zenitani & Hoshino 2007). 

How can such a system be simulated on a computer? Due to a relatively limited dynamic 
range of simulations, the outer boundaries will tend to affect strongly the injection region, 
so that the rate of injected magnetic flux would tend to adjust to the expansion rate. As we 
argued, this is physically unrealistic. In simulations the (anomalous) resistivity should be 
adjusted so that all the injected flux (or at least most of it) is destroyed before it is reflected 
from the outer boundary and affects injection. 

Komissarov (2006) gives a numerical example of resistive force-free magnetosphere, 
which has many similarities to the proposed model. Due to strong resistivity in the equato- 
rial current, the force-free wind converges to the equator, so that all flux surfaces formally 
become closed. The observed behavior seem to be independent of the resistivity, suggest- 
ing that a small numerical resistivity is sufficient to relax to a particular solution. The 
electromagnetic energy flows into the current sheet and disappears inside of it. 
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The inertial and kinetic effects will modify the structure of the equatorial current sheet 
and polar current. Conversion of magnetic energy into heat and plasma pressure will slow 
down the influx toward the reconnection regions. On the other hand, increased pressure 
will lead to faster radial expansion in those regions and will reduce this effect. Simulations 
of Komissarov (2006), do show that some flux surfaces close outside the light cylinder in 
resistive MHD simulations. 

The original impetus for this work came from discussions with Roger Blandford. I 
would like to thank Elena Amato, Jon Arons, Oleg Kargaltsev, Yuri Lyubarskii, Christopher 
Reynolds and Mallory Roberts for comments and discussions. 
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A. Separable solution of Eq. (13) 

Eq. (13) also allows a separable solution with distributed currents and charge densities 

7,*oc(l-^)sin 2 (Al) 

The radial component of the Poynting flux is oc (r — R 5 /r 4 ) 2 cos 9 sin 2 9, with flow lines given 
by sm 2 2 g = (-R 5 / r o — ^o) /(R 5 / r ' 2 — r ' 3 )- All the flow lines in this solution end at the equator. 



